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1. Introduction 



Precise calculations of the QCD cross-sections for the LHC are based on 
the factorization procedure. It allows to divide the entire collision process 
into separate parts: (1) the non-perturbative initial distribution of par- 
tons in the incoming protons, to be taken from experiments in form of the 
Parton Density Functions (PDFs); (2) the actual Parton Showers (PSs), 
both in the initial and final states, described by the evolution equations re- 
sulting from the resummed perturbative calculations; (3) the hard scatter- 
ing calculated perturbatively to a fixed order and (4) the non-perturbative 
hadronization of partons, described by phenomenological models. There are 
a few approaches to the factorization, amongst which the collinear factor- 
ization, well founded in the field theory, is the most popular one. Let us 
recall a few milestones in its development. The first, LO, resummations, 
dated to the early 1970s, are due to Gribov, Lipatov, Altareli, Parisi and 
Dokshitzer [I]. The general theorem has been addressed in the papers of 
Ellis, Georgi, Machacek, Politzer and Ross [2] with the help of the axial 
gauge in which the appealing ladder-like structure of the singularities has 
been proven. The practical calculations of the NLO kernels were done by 
Floratos, Ross and Sachrajda pip] and then by Kalinowski, Konishi and 
others [5j[6], based on the operator product expansion as well as by Curci, 
Furmanski and Petronzio [7j[8] in diagrammatic approach in dimensionally 
regularized MS scheme. Further reformulations of the factorization scheme 



were due to Collins, Soper, Sterman [9], Bodwin 10 and others in 1980s. 
Finally, twenty years later, the NNLO kernels have been calculated by Moch, 
Vermaseren and Vogt |11[|12| . All the above papers concern the analytical 
results. The Monte Carlo (MC) implementations, i.e. the PS programs, used 
in the actual data analysis were developed in a different pace. The first im- 



plementations: PYTHIA [13] and HERWIG 14 , based on the (improved) 
LO calculations, were created in mid-1980s and systematically developed 
since then. The next step - the complete NLO-based simulations have not 
been constructed yet, despite the fact that the analytical results are known 



for more than 30 years! There are two half-way solutions: MC@NLO 15 
and POWHEG |16| which combine the NLO hard matrix element with the 
LO-type PS (LO evolution). So far a LO accuracy of the MC programs was 
satisfactory in the data analysis, as compared to the experimental precision 
of the hadronic experiments. However, with the start-up of the LHC the 
situation has changed. The generic expected experimental precision of LHC 
results is of the order of 1%, and to match this precision the MC simulations 
must be of the NLO type at least. 

The goal of the KRKMC project [17, 18 developed in Krakow is to fill 
in precisely this gap and to construct for the first time the complete NLO- 
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Fig. 1. The collinear factorization for the Drell-Yan type process 

level PS program for QCD. In the following, after a short introduction to the 
collinear factorization, we will briefly describe the four main ingredients of 
the KRKMC project: the extended, exclusive collinear factorization scheme, 
the new LO PS scheme, exclusive NLO corrections to the hard process and 
exclusive corrections to the ladder part. 

In the results presented here some simplifications are temporarily present: 
only non-singlet kernels are included in the evolution, only qq — > W/Z hard 
process is included and non-running as is used. 



2. Extended exclusive collinear factorization scheme 

The standard collinear factorization rearranges the squared matrix el- 
ement into a sum of two-particle-irreducible universal kernels Kq, and a 
process-dependent function Co, as shown in Fig. [T] for the Drell-Yan type 
process (use of the physical gauge is instrumental here!). The collinear sin- 
gularities are located in the Kq functions. At this point one introduces 
the projection operator F which separates the singular part P-fTo from the 
non-singular rest (1 — W)Kq, leading to the rearrangement 

\M\ 2 = Cv—^r = C ®{l + K + K®K + K®K®K +...), (1) 
1 - A 

C = C °- 1-(1-P).AV K = FK °- 1-(1- 1 P). Ko < 2 > 

where the reorganized kernel and the hard process, K and C, are to be 
calculated up to a requested accuracy: LO, NLO, etc. The dot (•) denotes 
four-dimensional integration, whereas the ® is a one-dimensional integral 
(convolution) over the longitudinal momentum (lightcone) variable x + . The 
transverse momentum degrees of freedom are integrated out. 

Why is the above scheme not good for precision MC simulation? The 
three main reasons are: (1) It does not conserve four momenta. In order 
to reduce • into <S>, the F operator allows for unphysical configurations 
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of momenta. (2) There are very strong over-subtractions (cancellations) 
present. Eq. is a geometrical series whereas we expect the final result 
to be of the exponential form. (3) The scheme is defined in dimensional 
regularization (pole-part extraction), whereas the MC simulation must be 
done in four dimensions. How can one cure these drawbacks? (1) The 
projection operator F has to be redefined. (2) The time-ordered exponential 
must follow directly from the construction of the modified factorization 
expression ([I]). (3) The geometrical regularization has to be introduced for 
the real emissions instead of the dimensional one. 

We define a new JP' operator and also introduce a finite-part operator 
B. To the lowest order we have Ib^o] =K - ^'^{Kq}, as expected. In 
higher orders B is defined recursively: Ib^Ko ■ K ] = K ■ K - P~^{ S2 #o}- 
P^K^-^'^Kq- t S2 [K ]}- fe^Ko]- P^{K }, and so on. 

1p ^ sets the upper limit \i on the phase space for all real partons towards 
the hadron using the kinematical variable s(ki, .., k n ) < [i where s stands for 
virtuality, maximal rapidity, maximal k T , etc. The notation { Si j4} defines 

a = Si , e.g. Ip'^A} ■ Ip'^A} ■ 1p' S2 { s 'A} means ^> S3 > S2>Sl instead of 
^{-s A} ■ Ip^A} ■ Ip^A} corresponding to ^ >S3 9^ >S2 6^ >S1 (CFP- 
like). Also, P~^(j4) extracts a singular part from an integrand A (not from 
the integral / A like CFP!). 

An exclusive parton density function is then defined as the integrand in: 



D(jj) = exp TO ^ !^K - B" 5 



1-Kn 



, d (t D(n,x) = P®D(n)(x), 

(3) 



and at the inclusive level it fulfills the standard DGLAP equation. 



3. LO Parton Shower revisited 

Having established the general factorization framework, we now turn to 
the LO MC. It is the basis of the project. All the NLO effects will be added 
on top of it by means of the rejection techniques. The strategy is standard 
- one combines twice the formula ^ (a Drell-Yan case): 

oo oo 

^rrfrg)) = ^ ^ {^(Px^r^K^r]}^ (4) 

rti = l ri2=l 

and parametrizes it in terms of the Sudakov variables (ati,/3i for emitted 
partons and X{ for virtual, "ladder", ones) that can be directly generated. 
There are two technical problems to be resolved: (1) the constraint on the 
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value of the final xp/B = Y\z F , B (necessary for resonant processes) and 
(2) complete coverage of the phase space without any gaps. The standard 
solution of the problem (1) is the "backward evolution" 13 , which uses a 
pretabulated grid of PDFs. In order to avoid this complication we proposed 
a different algorithm which imposes the constraint on top of the normal 



"forward" evolution, see [19 ,20 for details. The problem (2) is solved by 
a kinematical mapping of the original phase space to the "t ang ent space" 
(ki,a{,Pi —7- ki,6(i,(3i). The mapping we proposed in Ref. tell is a plain 
rescaling, although defined in a recursive way: 

kir^Xikn, \ = — r J ' * = l > 2,...,ni + n 2 . (5) 

The rescaling factors Aj are chosen such that 



Si = sxi = s Y[ zfttj Y[ = ( P ~ kw ^ 2 = ( P ~ A J^j) 2 - ( 6 ) 

ttjGF TT 3 eB j=l 3=1 

This mapping has three important features: (1) it preserves angles, i.e. 
rapidity ordering, and the upper limit of the phase-space (integration) in 
rapidity; (2) it preserves soft factors (da ja . . . ), i.e. the original behavior 
in the soft limit; (3) it covers the phase-space completely, without any gaps. 
The complete LO MC algorithm looks as follows: 

1. the variables zp and zb are generated by the FOAM MC Sampler 



22,23 



2. the four-momenta are generated separately in the F and B parts of 
the phase space with the constraints Ylj&F ®j = 1 — -£f and J2jeB Pj = 
1 ~ zb, 

3. the double-ordering permutation ir is established, 

4. the rescaling parameter Ai is calculated; k n = X±k m is set, such that 
(P - k ni ) 2 = sx u 

5. the parameter A2 is calculated and k W2 = ^k^ is set, such that (P — 
k 1Vl — k^) 2 = SX2 = sz ni z n2 and so on, 

6. in the rest frame of P = P — k nj four- momenta and q% are 
generated according to the Born angular distribution. 



6 



USTRON2011 PRINTED ON NOVEMBER 24, 2011 



Exact analytical integration of the LO MC distributions of eq. Q over the 
multigluon phase space is possible (we use the rapidity ordering with 5 be- 
ing the rapidity of the produced boson, or equivalently the splitting point 
between the forward and backward hemispheres): 

/■I 

o-(C^r^r^) = / dx F dx B D F (E,x F ) D B (E,x B ) a B (sx F x B ), 
Jo 

with two PDFs obeying the DGLAP non-singlet LO evolution equation: 
~D F (E,x) = [P®D F (E))(x). 



4. Exclusive NLO corrections to the hard process 



Having described the "underlying" LO MC, we proceed now with the 
NLO corrections. We begin with the hard process part. The NLO correction 
consists of a one-parton real emission and the matching first order virtual 
correction. It is included by means of a MC weight. In the case of the 
Drell-Yan process (Fig. g the weight reads: 



W NLO. 
W MC - 



-, . A , sr^ (3i(s,p F ,p B ;aj,z Fj ) 

1 + <-±s+V + > — 

f£.P(z Fj )dv B (i,6)/dSl 



f3i(s,p F ,p B ;aj,z Bj ) 
'■^P(z Bj )da B (s,6)/dn' 



The IR/colinear-finite real emission part /3\ (numerators) is a fully-different- 
ial real emission matrix element minus two counter-terms (the counter-terms 
are LO collinear distributions generalized to off-collinear regions): 



Pi(PF,PB\qi,q2,k) 



l + (l-a-(3) 2 da B ^ g 

9a> ? 2 dn n {s ' 9) 



[l-a) 2 da B {l-P) 2 da B 

[s,e F i) + -j^-is^m) 



dCl q 2 dflq 

, l + (l-a-P) 2 da B ^ s 

a< ? 2 dn„ {s > 0) - 



The denominators are just the LO "underlying" differential distributions: 
numerators of the DGLAP LO kernels P{z F j/ B j) times the Born cross sec- 
tion da B (s,6)/dtt. The sums go over all emitted partons in both the F and 
B hemispheres. This way all partons contribute to the hard scattering and 
there is no problem of defining the "last" , "hardest" , etc. one. 

Similarly to the real emission part, the virtual+soft correction is defined 
as a difference between the complete inclusive result and two inclusive MC- 
type counterterms 



A 



v+s 



Dffi{z) - 2C™ c (z) 



C F a s 



<5-part 



7T 



-TT 
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D%${z) can be taken from the literature (eq. (89) in Ref. [24]) or re- 
calculated from the Feynman graphs [25]. Details, in particular on the 
construction of the counterterms, can be found in [21] . Note that the Ay+s 
is kinematics independent (it is a constant). 

As in the LO case, the exact analytical integration of the NLO MC 
distributions over the multigluon phase space is possible, and leads to: 

i 

a^Q^TpTB) = J d%F oIxb dz Dp{E,xp) Db{£,xb) ctb(sz%f%b) 
o 

x {5 z=1 (l + A s+ v) + C^ c (z) 



where 

« = --(i-4 (7) 



7T 



The above inclusive NLO correction to the hard scattering in the MC scheme 
differs from the MS correction (eq. (90) in [24]) 



C£ 8 (z) = C ^f^ [41n(l - z) - 2 In,]. (8) 

As we see, the MC correction in eq. Q is a simple, regular polynomial, 
whereas the MS one in eq. ^ contains the singular logarithmic terms ln(l — 
z) /(1—z). These terms originate from the "mistreatment" of the phase space 
done by the kx ordering implied by MS, see also [24] . 

To summarize, the new NLO MC factorization scheme has the following 
interesting features: (1) The NLO corrections are added on top of the LO 
MC with a simple, positive weight. (2) There is no need to correct for the 
difference in the collinear counter-terms of the PSMC and MS schemes pro- 
vided PDFs are in the PSMC scheme. (3) The virtual+soft corrections are 
completely kinematics independent - all the complicated (f £ c=t contributions 
of the MC@NLO scheme are absent. (4) There is a built-in resummation of 
the ln }^. x ^ terms (demonstrated for n = 1). 



5. Exclusive NLO corrections to the ladders 

The NLO corrections to the multi-emission ladder part have not been 
included in any of the existing PSs yet, except for a partial proposal of 
Ref. 1 26]. The scheme presented here is the first complete solution for the 
non-singlet evolution. Our scheme is based on the reweighting technique. 
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Once again, we begin with the LO "underlying" differential distribution for 
the single ladder: 



D LO (x,Q)=J2 



n=0 



oo n ,3 , 

^ E 1 1 " 

n=0i=l ^* 



^Q>a i >a i - 1 Pl (ki), 



2C 2 F a s 1 1 + z 2 



7T 



T2 



The NLO corrected distribution looks as follows: 



D 



NLO, 



x,Q) 



E 

n=0 



n pi— 1 

E E 

pi=l 3i=l 



oo / n 



^+e n 



n=l v j=l 



n pi — 1 Pi-l P2 — 1 

+E E E E 



Pl=lP2 = l Jl =1 J2 =1 



A- 



+ 



1+ 



Q>a i >a i _ 1 

n p— 1 n Pi— 1 pi— 1 P2— 1 



p=l 3=1 



Pl=lp2=l Jl = 1 J2 =1 

il#P2 J2#P1J2 



The above formula might look complicated, but in fact its structure is sim- 
ple. The LO ladder is multiplied by appropriate non-singular NLO weights: 
in the first term by the virtual weights /3q (in the picture dots are re- 
placed by squares); in the second term one NLO real correction, W(k p , kj), 
is added with all possible choices of the first and second momentum (two 
sums are visible); the third contribution is identical to the second one but 
with two NLO corrections and so on. The actual definitions of the weights 
are intuitively obvious - they are ratios of the NLO to LO distributions 
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(with counterterms if needed) and can be graphically represented as: 



i + 2»(a; 1 . 



ISR' 



W(k 2 ,h 



m 2 m\ 



i. 



We have recalculated the real emission weight W(k2, k\) in the new factor- 
ization scheme 27,28 and we are in the process of recalculating the virtual 
weight |25| . We have tested numerically the above algorithm by compar- 
ing the inclusive distribution D NLO (x,Q) from our exclusive MC with the 
standard inclusive result. We obtained three-digit agreement limited by the 
statistics, see [T7J[l8] for details. 

In the above we have ignored the contribution to W(k2,ki) from the 
gluon-pair production. This contribution has additional singularity in the 
limit of a vanishing mass of the emitted gluon pair. In the inclusive ap- 
proach this singularity is cancelled by the corresponding virtual one. In 
the MC exclusive simulation this is not possible and the singularity has to 
be included into the "underlying" distribution, otherwise it would ruin the 
convergence of the algorithm. Fortunately, as it is well known, this singu- 
larity is a part of the LO corrections to the final state and can be included 
into the algorithm at the expense of a few additional combinatorial sums 
over the final-state emissions. Graphically it looks as follows: 



SR ^FSR 



oo m 

£ £ 

n,m=0 r=l 



■fr 4- --4- ■».■!]-■ 



n-l 
'n-2 



where the Sudakov form-factor S 

2 



, R is subtracted in the virtual part: 

2 



{l + 2K(A ISR + V FSR 



'FSR 



and the FSR real counterterm is subtracted together with the ISR one: 

Both the virtual and real corrections (weights) are now regular. All singu- 
larities have cancelled separately in each weight, allowing for construction 
of the efficient MC algorithm. The exact analytical integration is possible 
in this case and we have used this result to perform numerical tests of the 
MC program with the percent-level precision. 
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6. Summary 

We have briefly reviewed the most important results of the first ever, 
complete NLO, fully exclusive, PS project under development in Krakow: 
(1) Extension of the collinear factorization, better suited for the MC im- 
plementation, is defined. (2) LO PSMC is (re-)constructed from scratch, in 
a way compatible with the new factorization scheme. (3) The NLO parts 
of the hard process and the evolution kernels are recalculated in the new 
scheme (so far the non-singlet NLO exclusive kernels only). (4) The dif- 
ferences between the new MC and standard MS schemes are understood, 
keeping the universality (process independence) in mind. (5) The proposed 
solution has advantages as compared to the other techniques of adding the 
NLO corrections to the hard process and it is completely new for the ladder 
parts. (6) Implementation in the MC has been tested at the prototype level 
with the relative precision of 10 -3 . 
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